Tuning the bond order wave (BOW) phase of half-filled extended Hubbard models 
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Theoretical and computational studies of the quantum phase diagram of the one-dimensional 
half-filled extended Hubbard model (EHM) indicate a narrow bond order wave (BOW) phase with 
finite magnetic gap E m for on-site repulsion U < U* , the critical point, and nearest neighbor 
interaction V c w U/2 near the boundary of the charge density wave (CDW) phase. Potentials with 
more extended interactions that retain the EHM symmetry are shown to have a less cooperative 
CDW transition with higher U* and wider BOW phase. Density matrix renormalization group 
(DMRG) is used to obtain E m directly as the singlet-triplet gap, with finite E m marking the BOW 
boundary V S (U). The BOW/CDW boundary V C (U) is obtained from exact finite-size calculations 
that are consistent with previous EHM determinations. The kinetic energy or bond order provides 
a convenient new estimate of U* based on a metallic point at V C (U) for U < U* . Tuning the BOW 
phase of half-filled Hubbard models with different intersite potentials indicates a ground state with 
large charge fluctuations and magnetic frustration. The possibility of physical realizations of a BOW 
phase is raised for Coulomb interactions. 

PACS 71.10.Fd, 71.30.+h, 71.45.Lr 
PACS numbers: 



INTRODUCTION 



The quantum phase diagram of the half-filled extended 
Hubbard model (EHM) in one dimension (ID) illustrates 
competition among on-site repulsion U > 0, nearest- 
neighbor interaction V and electron transfer t. Large U 
gives a ground state (gs) with one electron per site while 
large V > leads alternately to empty and doubly oc- 
cupied sites. The gs phase diagram and its critical point 
U* have evolved since Hirsch's original study [lj of the 
boundary between a spin density wave (SDW) at large 
U and a charge density wave (CDW) at large V. Naka- 
mura [2fl first proposed a bond order wave (BOW) phase 
between the SDW and CDW up to a critical U = U*. 
The schematic quantum phase diagram in FigJT] follows 
Sengupta et al. [3| and recent works by Zhang [J] and 
by Glocke et al. 0. The CDW boundary that we de- 
note as V C (U) has readily been found numerically. The 
value of U* and the SDW boundary at V S (U) < V c 
are more challenging. The BOW phase is much nar- 
rower than sketched and lies between 0, [f| V s and V c for 
U < U* « It. It is insulating and has a finite magnetic 
gap, E m , that is expected to be very small because V S (U) 
is a Kosterlitz-Thouless transition [2j . Finite E m implies 
electronic dimerization leading to long-range order at 
K. 

Multiple theoretical methods have been applied to the 
quantum phase diagram at K. They include field- 
theoretical analysis of continuum models 0, @, weak- 
coupling expansions Q , Monte Carlo simulations 0, [1] , 
density matrix normalization group (DMRG) methods 
[H B @] and the transfer matrix renormalization group 
(TMRG) @, in addition to exact diagonalization [10( 
of finite systems and perturbation expansions about the 



weak and strong coupling limits. The BOW phase for 
weak coupling (small U, V) has recently been demon- 
strated Q for the EHM by functional renormalization 
group. 

In this paper, we study the K phase diagram of 
half-filled extended ID Hubbard models with intersite 
interactions V m that are not restricted to V\ = V. Any 
spin-independent V m retains [llj the translational, spin, 
and electron-hole symmetry of the EHM. Suitable po- 
tentials have electrostatic energy — Vcym/2 per site, with 
Madelung constants <xm = 2 for EHM, 21n2 for point 
charges, etc. The CDW transition at large V becomes 
less cooperative for 1 < a a/ < 2, thereby extending the 
BOW phase to larger critical U* and increasing the width 
V C (U) — V S (U). We introduce a convenient new way to 
find U* and V C {U) based on the kinetic energy or bond 
order p(U, V), which is maximized at a metallic point 
where the large U and V are balanced. The threshold 
of the BOW phase for U < U* is obtained by DMRG 
calculation of the magnetic gap E m . Within the DMRG 
accuracy, finite E m at V S (U) below V C (U) sets a lower 
bound on the onset of the BOW phase. 

Competing interactions are widely invoked for materi- 
als whose gs depends sensitively on small changes [l2j]. 
The EHM literature deals with many-body techniques 
and ideas without regard to possible physical realiza- 
tion. Coulomb interaction between point charges is more 
physical. As in metal-insulator or Mott transitions, how- 
ever, the assumption of purely electronic transitions is 
inherently an approximation. The SDW is unstable to 
a Peierls transition leading to dimerization for any har- 
monic lattice. Lattice dimerization corresponds to well- 
understood BOW phases that are realized in quasi-lD 
systems such as organic ion-radical salts [l3| and charge- 
transfer (CT) salts [H| or conjugated polymers flol fig. 
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FIG. 1: Schematic quantum phase diagram of half-filled ex- 
tended Hubbard models with on-site U > 0, nearest-neighbor 
V and ID potential with Madelung constant au- The 
SDW/CDW transition line V C {U) is first order for U > U* , 
the critical point, and second order for U < U* . The BOW 
phase between V s (U) and V C (U) is a Mott insulator with a 
metallic point at V c and a finite magnetic gap, E m , that opens 
at Vs. 



Such systems have intermediate correlation (U—V) > 2t. 
The present study is restricted to a rigid ID lattice with 
purely electronic instabilities. Possible physical realiza- 
tions are mentioned in the Discussion. 

Previous EHM studies [E S, S i, 1, S, 0, S & M 
have focused on gs correlation functions for charge, spin 
or bond order as well as various susceptibilities. Our 
DMRG calculations target instead E m directly as the 
singlet-triplet gap, Est, between the lowest triplet state 
and singlet gs. We also present exact finite-size results 
that have been used to model the neutral-ionic transi- 
tion of CT salts [ItJ • Contributions due to t are particu- 
larly important at V C (U). The line V C (U) up to U* that 
marks a continuous transition is a metallic point with 
special properties according to the Berry-phase formu- 
lation of polarization [ljl The motivations of the 
present work are to increase the range and width of the 
BOW phase by changing the intersite potential. In ad- 
dition to the computational advantages of a wide BOW 
phase, control of the width is a step towards understand- 
ing a poorly characterized gs. On the SDW side, we refer 
to previous mappings [20j, |21| of Hubbard models onto ef- 
fective spin Hamiltonians with frustration [22] . Different 
perspectives on the K phase diagram of EHMs are con- 
sistent with other EHM results that also serve as checks 
for long range potentials. 

The paper is organized as follows. Extended Hubbard 
models with spin-independent intersite potentials are de- 
fined in Section II along with bond orders, phase bound- 
aries and the metallic point. DMRG calculations for the 
magnetic gap E m and exact finite-size results for V C (U) 
and U* are presented in Section III and compared to 
prior EHM results. The range and width of the BOW in- 



crease for less cooperative potentials. Section IV briefly 
discusses possible realizations of BOW systems and inter- 
pretations in terms of charge fluctuations and magnetic 
frustration. 



II. PHASE BOUNDARIES OF HUBBARD 
MODELS WITH EXTENDED INTERACTIONS 

We consider a half-filled ID Hubbard model with pe- 
riodic boundary conditions (PBC), interaction potential 
V m between m th neighbors, on-site repulsion U > and 
nearest-neighbor transfer t, 



p. a 

p rn>0 



i) 
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(i) 



The number operators h p have eigenvalues n p = 0, 1 or 2. 
As written, the interaction energy is zero when n p = 1 
at all sites, at density p = 1. The Hellmann-Feynman 
theorem gives the weight or density of doubly occupied 
sites in the gs as 



P2 (U,V) = (2N)- 1 J2(n P (n P ~l)) = 



1 dE 
N dU 



(2) 



where Eq is the exact gs energy. Stoichiometry relates the 
densities of cationic sites (holes) with n p — and anionic 
sites (electrons) with n p = 2 as po = p-i = (1 — p)/2. The 
choice of V rn is open, subject to the constraint of spin- 
independent interactions with V = V% . We note that 



V m (a) = V 



+ 

(a + m) 



(3) 



corresponds to the EHM in the limit a — * — 1, to a point 
charge model (PCM) at a = 0, and slower decrease for 
a > that mimic molecular sites in a delocalized charge 
model (DCM). The electrostatic energy per two sites of 
the CDW with holes on one sublattice and electrons on 
the other is 



E M (a) = 2 (~l) m V m (a) = -Va M (a) 



(4) 



m>0 



The Madelung constant olm defined in Eq. |4] is easily 
evaluated analytically for integer values of a. 

In the strong-coupling or localized limit (t = 0), the 
SDW/CDW boundary is simply U = VctM, the dashed 
line in Fig. [TJ The gs for large U has n p = 1 at all sites 
and two-fold spin degeneracy at each site, while the gs 
at large V is a doubly degenerate CDW. The curvature 
23] of V m (a) ensures that the gs at t = is either the 
SDW or the CDW. In this limit, we have dV/dU -> l/a M 
and the density p 2 jumps from to 1/2 at V c — V jau on 
increasing V at constant U . For finite t, the discontinuity 
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of p2 at V C (U) decreases and vanishes at V* = V C (U*) 
when the transition becomes continuous. 

Translational symmetry leads to uniform bond order 
p(U,V), proportional to the gs expectation value of the 
kinetic energy part of Eq. [T] 



P(U, V) = — ^2 (( & Ua P +ia + h.c.)) = - 

pa 



1 dE 
2N dt 



(5) 



U = V = is a metallic point with po = 2/ir for a half- 
filled ID band. The bond order of extended Hubbard 
models has a maximum at V C (U). A first-order transi- 
tion for U > U* is indicated by discontinuous p(U, V) at 
V C (U), while a second-order transition for U < U* has 
continuous p(U, V) that develops a kink at p(U*, V*). 

The modern theory of polarization in solids grew out 
of the necessity of incorporating PBC for practical super- 
cell calculations [HI]. Subsequent generalizations of the 
theory and its relation to Berry phases have diverse ap- 
plications, including how to distinguish between metals 
and insulators [l9[ . The theory is applicable to correlated 
models with neutral-ionic transitions [24] or to models in 
Eq. [1] with any potential V m . For EHMs, we require 
the exact gs of ID supercells with N sites and PBC to 
compute the expectation value [2J] 



Z N (U,V) = (exp(^)) 



(6) 



M is the conventional dipole operator for unit charge and 
unit spacing in ID 



M = ^p(n p - 1). 



(7) 



Extended models in Eq. [T] have inversion symmetry at 
each site and at the center of each bond. The correspond- 
ing symmetries for finite N are reflections through sites 
or through bonds. Either symmetry ensures that Z/v is 
real in general for EHMs and reduces Zjq to the twist 
operator [l8[ with cos(2ttM /N) instead of the exponen- 
tial in Eq. [6] Quite generally, we have Zn ~ 1 when 
the gs consists largely of sites with n p — 1 in the SDW 
phase and Z^ w — 1 when the gs is a CDW with n p — 2 
on one sublattice, n p = on the other. Hence the sign 
of Zyv changes as V increased at constant U. The point 
Zn(U,V c ) = for U < U* corresponds to a metal [HI 
that separates two insulating phases. 



III. MAGNETIC GAP AND METALLIC POINT 

In this section we study the gs of H in Eq. Q]for poten- 
tials V m (a) using density matrix renormalization group 
(DMRG) and exact finite-size calculations. Since H con- 
serves total spin, E m is the singlet-triplet gap between 
the singlet gs and the lowest triplet state. We take t = 1 
as the unit of energy and find E m as the difference be- 
tween the lowest energy states with total Mg = 1 and 0. 



DMRG with PBC is found to be very accurate when the 
ring is expanded from the middle by two sites at each 
step. Earlier DMRG introduced one site in the middle 
and another at the end of the block spins [25| . This had 
the disadvantage of adding a bond between a new site and 
the first site, whose operators have already been renor- 
malized many times. The accurate DMRG-PBC proce- 
dure treats the ring as two chains that are joined at the 
ends. 
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FIG. 2: Infinite DMRG procedure for a system with periodic 
boundary condition. Sites on left and right blocks are num- 
bered as unprimed and primed integers respectively. Old sites 
are shown as filled circles and sites added at the DMRG steps 
are represented by crosses. 

At each DMRG step, two sites are added, alternately, 
in the middle of one chain. The schematic diagrams are 
shown (see Fig. [5]). The new transfer terms at 2n system 
size are (n, n'), [(n — 1), (n — I)'], [(n, (n — 2)'] 
and [(n — 2), n'\. The electron-electron interactions are 
diagonal and hence are known to be accurate even when 
interacting sites are separated by several renormalization 
steps. The procedure introduces explicit transfers be- 
tween sites whose operators have been renormalized only 
twice. We have carried out finite DMRG calculation for 
every An system size for increased accuracy. 

To calculate E m , we retain to = 150 dominant density 
matrix eigenvectors for states in Eq. [T]with either total 
M s = or 1. Each DMRG step gives E m (N) with N 
increasing by two. We have performed DMRG calcula- 
tions with finite ring size N = An and find that one finite 
DMRG sweep is sufficient for good convergence of the en- 
ergy. Fig. [3]shows representative E m (N) vs. 1/N results 
up to N > 50 for the EHM (a = 1) and PCM (a = 0) 
at U = 4 for several values of V. Good 1/N dependen- 
cies are found and yield the extrapolated E m in Fig. [3] 
Other U and potentials V m (a) show similar 1/N behav- 
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ior for small E m when V < V C (U), while V > V c results 
in large E m « 1 that increase rather than decrease with 
TV. We estimate the accuracy of extrapolated E m 's by 
comparison with E m = in the Hubbard model (V = 0) 
or in the SDW phase (V < V S (U)). The extrapolated 
E m /t are « 0.005. Accordingly we have assumed that 
the spin gap is finite when E m /t > 0.01. We verified that 
E m does not depend significantly on the choice of m for 
m > 130. Much larger m leads to computational diffi- 
culties since the sparseness of the Hamiltonian matrix is 
significantly reduced. The 1/N dependence of E m {N) is 
monitored to find N after which the extrapolated value 
changed by less than 0.001, and another iteration was 
then performed. 



LU 



1.0 r 

0.8- 

0.6 

0.4 

0.2 

0.0- 
1.: 



U=4 




0.16 
0.14 
0.12 
0.1 

E 

W0.08 
0.06 
0.04 
0.02 
0, 



EHM U=4 



# V=1.7,E =0.004 

' m 

a V=1.8,E =0.004 

m 

♦ V=1.9, E =0.023 

' m 

■ V=2.0, E =0.029 




o V=1.7,E m =0.001 

a V=1.8,E™=0.003 

O V=1.9,E™=0.005 

□ V=2.8, E m =0.048 







0.005 0.01 0.015 0.02 0.025 0.03 

1/N 



FIG. 3: Singlet-triplet gap, E m , of Eq. [Tjwith TV sites and 
periodic boundary conditions vs. 1/TV for EHM with a = — 1 
in Eq. |3]and PCM with a = 0. Extrapolation is between TV = 
28 and TV > 50. The insets specify U , V and the extrapolated 
|£ m |<0.01. 

To show the threshold V S (U) for opening a magnetic 
gap, we plot E m in Fig. H at U = A for EHM and PCM 
(a = in Eq. [3J. Dashed lines V s at E m = 0.01 indi- 
cate the SDW/BOW boundary within our DMRG accu- 
racy. The magnetic gap opens more slowly for PCM than 
for EHM and more slowly still for delocalized charges 
(DCM, a = 1, not shown). Dashed lines at V c indicate 
the BOW/CDW boundary that is found below. Similar 
DMRG evaluation of E m — 0.01 for other U and V m (a) 
combinations yield the V S (U) entries in Table 1. We find 
the EHM gap to open at V s = 1.86 for [7 = 4. In this 
case, direct evaluation of E m is close to early estimates 
0, 0] of the opening of the magnetic gap, while recent 
calculations 0, Q give V s « 2.02. The DMRG result in 
ref. 4 is based on a broad peak at V s that would pre- 
sumably sharpen in larger systems. The TMRG result in 
ref. 5 is a thermodynamic method at finite temperature 
that nevertheless gives estimates for K properties on 
extrapolation. We find a substantial E m — 0.036 for the 
EHM at U = 4 and V = 2. 

The CDW boundary at V C (U) has been found by di- 
verse methods, as tabulated in ref. 8 for the EHM at 



FIG. 4: Extrapolated singlet-triplet gaps, E m , vs. V for EHM 
and PCM at U = 4. The vertical dashed lines V s marks 
E m — 0.01, the DMRG uncertainty (see text). The vertical 
dashed lines V c mark the metallic point 14(4) discussed in the 
text. 



[7 = 4; the range is 2.10 < 14(4) < 2.16. Here we eval- 
uate V C (U) for finite systems of N sites and PBC in Eq. 
[TJ Exact correlated states of Hubbard- type models, cur- 
rently up to TV = 16, are found using a many-electron ba- 
sis of valence bond diagrams [2(| . As discussed in connec- 
tion with the neutral-ionic transition in donor-acceptor 
stacks the gs has a symmetry crossover at V C (U, N) 
for TV = An that depends weakly (» 1/TV 2 ) on TV. The 
bottom panel of Fig. [5] shows the U = 4 crossover of the 
EHM from a singlet gs with even electron-hole symmetry 
that transforms as k = tt for V < V c (4, 16) = 2.10 and 
as k = for V > V c . Moreover, we find pi(U, V) = 1/4 
and Zn(U, V) = for the k = singlet at almost exactly 
V = V C (U, TV) for U <U*. The EHM crossovers at U = 4 
for TV = 8, 12, and 16 extrapolate to 14(4) = 2.16 ±0.01, 
the entry in Table 1. Tighter extrapolation is possible 
using TV = An + 2 rings with antiperiodic boundary con- 
ditions [13, but this was not pursued. Finite-size correc- 
tions to 14 are even smaller for EHM at U — 6 or 10. 
They are somewhat larger for PCM, whose V C (U) entries 
in Table 1 have estimated uncertainties of ±0.02, and 
larger still for DCM with V C (U) uncertainties of ±0.03. 

The boundaries of the BOW phase for several poten- 
tials are listed in Table 1. The SDW/BOW threshold 
V S (U) corresponds to E m = 0.01. The BOW/CDW 
boundary V C (U) is the extrapolated gs crossover, the 
metallic point for U < U* . The U = 4, 6 and 10 en- 
tries clearly show increasing 14 — 14 as a m decreases 
and the CDW transition becomes less cooperative. The 
demonstration of a BOW phase is then less demand- 
ingthan for the EHM, where conflicting U = A results 
US Hi were debated until recently. In particular, sub- 
stantial E m /t > 0.2 are achieved for the PCM and DCM 
at V < 14(4). Such gaps are robust numerically and 
easily exceed ksT in organic 7r-radical stacks [l3[ with 
t w 0.1 — 0.3 eV. It has been fully appreciated that an 
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FIG. 5: Lowest energy singlets, both with e-h symmetry +1, 
for EHM with U = 4 and 16 sites. Open squares, B u (k — n); 
stars, A g (k — 0). Panels (a) and (b) are Z in Eq. 6 and pi in 
Eq. 2. Panel (c) shows the gs crossover, with E = for B u 
and AE for A a . 



exponentially small E m poses numerical difficulties. The 
present results are upper bounds for V s at finite E m . The 
BOW phase narrows at U — 6 and is absent for the EHM 
at U — 10. Our interpretation, continued in the Discus- 
sion, is that a BOW gs requires large t in Eq. [T] and 
strong charge fluctuations. 

TABLE I: BOW phase between V c — V s for three potentials, 
EHM, PCM and DCM, with a = -1, 0, and 1 in Eq. O at 
U — 4,6 and 10. The threshold V s (U) has magnetic gap E m — 
0.01, while V C (U) is the CDW boundary. The critical point 
U* has V s = V c and no BOW phase. The Madelung constant, 
olm in Eq. 3J decreases for less cooperative transitions 



Models 


EHM (a = -1) 

V 3 V c 


PCM (a = 0) 

V s Vc 


DCM(a = 1) 

Vs V c 


£1 = 4 


1.86 2.16 


2.50 3.02 


2.70 3.35 


U = 6 


3.06 3.10 


4.26 4.45 


4.63 5.03 


U = 10 


5.11 5.12 


7.33 7.42 


8.23 8.35 


u* 


6.7 ±0.2 


10.6 ±0.3 


13.3 ±0.5 


CtM 


2 


21n2 


4(1 - ln2) 



We turn next to the critical point U* , whose evaluation 
has been challenging. The most recent EHM values are 
U* = 6.7 ± 0.2 (ref. 5) or 7.2 (ref. 4), almost twice the 
first estimates from Monte Carlo simulations. We obtain 
U* from the kinetic energy at V C (U) as follows. The bond 
order p(U, V) in Eq. O is a continuous function of V for 
U < U* with a maximum at the metallic point V C {U). 
The peak narrows with increasing U. A kink develops 
at p(U*,V*) and p(U,V) is discontinuous for U > U*, 
where the transition is first order and there is neither a 
BOW phase nor a metallic point. Large U generates a 
discontinuity mp(U, V) at V C (U), with less kinetic energy 
on the CDW side. Second-order perturbation theory for 
the energy and Eq. [S] provide a simple expression for the 



bond order, 

It 

*W = vj aH -i) + <y-v M 2a u -ir v>v " 

Although U w 10 is not in the big-L/ limit, Eq. [H] 
rationalizes direct evaluation of p(U, V) and the onset of 
a discontinuity at U = V c (Xm- 

Figures [5] and [7] show p(U, V) of N = 16 rings for the 
EHM and PCM, respectively, for the indicated values of 
U . The band limit (U = V = 0) of the extended system 
is p = 2/7r; the analytical result for N = 16 is 1.3% less 

, . 1 ± V2 ± 2(coso!) ± sin0) 

po(16) = -^-^ — = 0-62842 (9) 

8 

where <fr = 2n/16. The band limit is the dashed horizon- 
tal line when plotted against (V — V c )/V c . In interacting 
models of Eq. 1, p(U, V) has a broad maximum at V c 
that is less than po and that sharpens with increasing U 
for either potential. The kinetic energy at the metallic 
point decreases slowly. Careful examination of the U = 6 
curve in Fig. [7]shows that for N — 16, p(6, V c ) is slightly 
larger for the gs on the CDW side. The [7 = 4 and 
6 curves in Fig. [S] show the same effect for the EHM. 
We used additional values of U as well as TV = 12 and 
16 results to estimate U* as the onset of discontinuous 
p(U, V). The U* entries in Table 1 increase as expected 
with decreasing aju- U* = 6.7 for EHM agrees with the 
two recent calculations [4, 5]. The dashed line in Fig. [6] 
for EHM is p^ in Eq. E|with V c = 5.12, scaled by a fac- 
tor of 0.80; the dashed line in Fig. 0for PCM is p {2) with 
V c = 11.65 and a scale factor of 0.75. The perturbation 
result captures most of the p(U, V) discontinuity. 
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FIG. 6: Bond order, p(U,V) in Eq. 5, vs. (V - V c )/V c for 
EHM with U = 4, 6, 7 and 10. The dotted line is the band 
limit, po = 2/7r; the dashed line is the strong-coupling limit, 
Eq. E]with V c = 5.1, t = 1, scaled by 0.80 
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Band limit (U=V=0) 
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FIG. 7: Bond order p(U, V) vs. (V - V c )/V c for PCM with 
U = 4, 6, 10 and 16. The dotted line is the band limit, 
po = 2/7r; the dashed line is the strong-coupling limit, Eq. [8] 
with V c = 11.65 and t = 1, scaled by 0.75 



The critical point Lf*,V* = K;(J7*) in Fig. Q] marks 
the BOW phase with a metallic point. The first-order 
transition for U > U* is between two insulating phases. 
One has pi < 1/4 and Z > 0, the other has P2 > 1/4 
and Z < 0, and there is no metallic point. In the limit 
of large U = aj\iV c (or small t), we have p ss according 
to Eq. [H] The behavior of p(C7, V) and its weak size 
dependence provides a convenient new way to find U*. 
The reason is that p(U, V) varies slowly with V near the 
maximum at V C (U), whereas p2{U,V) or Z(U,V) vary 
the most rapidly there. It is then numerically easier to 
discern a discontinuity in p(U, V). 

Half-filled Hubbard models with spin-independent po- 
tentials have e-h symmetry. Their dipole allowed (one 
photon) and two photon excitations are consequently to 
different manifolds of states, as discussed extensively for 
linear polyenes 27]. The lowest two-photon state, 2 1 A g , 
has the gs symmetry and can be viewed as two triplets. 
Then E m — on the SDW side implies that the two pho- 
ton gap also vanishes, while the one-photon excitation to 
l 1 i? u has finite energy. Conversely, finite E m in the BOW 
phase implies a finite two-photon gap. While DMRG for 
E(2 1 A g ) is less accurate at present than for E m , we find 
a two-photon gap of 0.138 for EHM at U = 4, V = 2.1, 
well within the BOW phase. Moreover, since the gap is 
much smaller than 2E m — 0.48, we conclude that the 
triplets form a bound state. A finite two-photon gap is 
another signature of a BOW phase. 



IV. DISCUSSION 

We introduced potentials V m (a) in Eq. [3] that retain 
the symmetry of the EHM in order to extend and widen 
the BOW phase sketched in Fig. [TJ DMRG calculations 
of the magnetic gap E m are a direct new approach to 
the SDW/BOW boundary. The kinetic energy or bond 



order p(U, V) at the metallic point turns out to be a con- 
venient new way to evaluate the critical point U*. Tun- 
ing the BOW phase by changing V m (a) also clarifies the 
electronic instability with increasing U < U* before the 
inevitable CDW transition at large V. The BOW phase 
is interesting but poorly characterized because it does 
not appear in the limit of either large U or large V. It 
is an intermediate phase with strong charge fluctuations, 
competition between U and V leading to magnetic frus- 
tration, and fairly large t. Translational and inversion 
symmetry are broken in the BOW phase, while trans- 
lational and electron-hole symmetry are broken in the 
CDW. Finite V is required in either case. 

The CDW shown in Fig. Q]is a good representation of 
the gs when 7>(7 and sites are alternately empty and 
doubly occupied. The SDW in Fig. Q] is a cartoon for 
U ^> V, however, since the gs of a Heisenberg antifer- 
romagnet (HAF) is a linear combination of many states 
with rip = 1 at all p. The Neel state shown has higher 
energy than the Kekule functions, \K1) or \K2) 1 in which 
adjacent spins are singlet paired, either 2p with 2p — 1 or 
2p with 2p + 1 for all p. \Kl) or \K2) is a linear com- 
bination of A-spin states, half up and half down, that 
includes the Neel state. The actual gs for U ^> V also 
has singlet pairing of sites that are not adjacent. 

We have deliberately omitted a BOW cartoon in Fig. 
B it is usually shown 0, |, 1, @ as (T |)(T I) to 
suggest broken symmetry for bond orders. Such a repre- 
sentation, suitable for a dimerized HAF, is quite mislead- 
ing for a BOW gs with almost 50% of empty and doubly 
occupied sites (Fig. [5]) and large bond order (Figs. [5] 
and [7]) for V slightly less than V C (U). Since t connects 
adjacent sites, large p(U, V) indicates strong mixing of a 
singlet pair at sites p, p+l and ionic singlets with n p — 2, 
n P +i = or rip = 0, n p+ i = 2. The BOW gs is a linear 
combination of the full Hubbard basis of A N states, aside 
from particle number or symmetry constraints. A prod- 
uct of dimer function (1,2)(3,4)... gives some insight into 
the BOW gs. The dimer function (1,2), an approxima- 
tion for the BOW gs, is 



(1,2) 



cos 9 
[ ~V2 

+ 7T 



0L4e - a \e a l a ) 
(<AAe + 4AeW) (10) 



where a) is a creation operator in Eq. [T] and 9 is a vari- 
ational parameter. The state in the first term is the co- 
valent (Heitler-London) singlet; that in the second is the 
ionic singlet. The bond order per site is sin 9 cos 9 w 1/2 
for equal admixture at 9 = tt/4. Such p(U, V) are found 
directly for U < U* . The product function indicates a 
BOW with large charge fluctuations 

Strong coupling gives a first-order SDW/CDW tran- 
sition at V c = U/atM f° r t — and no BOW phase for 
U > U* . Corrections for finite t can readily be found. 
At the boundary, the energy for an electron transfer in 
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either the SDW or CDW gs is 



AE = U - V c = U(l - 



a M (a) 



(11) 



We have AE > for cxm > 1 and can apply perturba- 
tion theory in t/ AE. Second-order {t/U) 2 corrections to 
E (U,V)/N on either the SDW or CDW side yield an 

approximate p 2 2 (U, V) according to Eq. [51 The jump at 
V r decreases and vanishes at 



2a M i\/2hi2 + 1 
{a M - 1) 



(12) 



ere) is a simple estimate for U* that shows how poten- 
tials with small c*m > 1 extend the range of continu- 
ous p 2 (U,V). We obtain U^/t = 6.18 for the EHM 
with a m — 2. Second order corrections yield a com- 
parably accurate estimate for a continuous neutral-ionic 
transition in CT salts [Hj]. The values of U {2) arc 11.1 
for PCM and 16.7 for DCM. Second order corrections 
yield V c (U,a) > U/ocm but overestimate the increase; 
the fourth order correction [2l[ for EHM reduces the in- 
crease. 

Half-filled Hubbard models have spin-charge separa- 
tion when U — V is large compared to t. The gs on the 
SDW side can be mapped into an effective spin Hamil- 
tonian for any potential V m . Up to t 4 , H e ff is a HAF 
with Ji, J2 > and frustration; a magnetic gap opens 
[22| at J2/J1 — 0.2411. The possible connection between 
a BOW phase and spin frustration has been pointed out 
previously [2lj . and it is a delicate matter. Seitz and 
Klein ||f obtained H eff up to {t/Uf for V = in Eq. 
[T] by systematically considering virtual transfers in the 
covalent sector with n p = 1 at all p; an even number of 
transfers is required to end up with n p — 1 at all p. Van 
Dongen [H[ obtained H eff for EHM up to [t/(U - V)} 4 . 
The t 2 term gives a HAF with nearest-neighbor exchange 
Ji = 4i 2 /([/ — V"), and J\ has contributions from all 
higher orders. The t order introduces a second-neighbor 
exchange J2, also AF, and all higher orders contribute to 
J 2 . Mapping into a HAF with J\ and J2 breaks down 
in the next order, even if such corrections are included 
in Ji and J2, because t 6 generates exchange interactions 
among four successive spins [2(| . At the t 4 level, the ratio 
Ji I J\ can readily be generalized to the potentials V m (a) 
in Eq. El We obtain 



J2/J1 



t 2 F(a,V) _ 

(U - V) 2 ~ 1 + 4J 2 / Ji 



F { V,a)-l= V - V ^ 



V 



U-V 2 (a) V + (a + 2)(U - V) 



(13) 



(14) 



The factor F(V, a) reduces to F = 1 at V = and to 
F(V, 1) = 1 + V/U for the EHM. Finite E m requires 
J2/J1 > 0.2411 in spin chains when the rhs of Eq. [TBI is 
0.1227. 



Since the ID Hubbard model is rigorously known [29| 
to have E m = for U > 0, F = 1 is insufficient to 
open a gap at U/t = 2.85 for V = in Eq. [TJ H eff 
with Ji and J2 clearly fails for (U — V) < 3i because 
it incorrectly predicts finite E m . The EHM has F — 
1.5 at V = U/2, when a magnetic gap opens at (U — 
V) ft = 3.5 according to Eq. [TJ and yields U* = 7.0 
in surprisingly good agreement with recent calculations 
d, [1] and Table 1. But the correction to J2/J1 goes as 
[t/(U — V)] 2 « 8% in the next order and -ff e // has an 
additional term. Potentials with a > 1 lead to F > 1.5 
at V = U/ocm and to slightly larger (U — V)/tPa 3.6 — 
3.7 in Eq. 1131 that, however, substantially overestimates 
U* in Table 1. Hence we regard the EHM result to be 
accidental. Growing magnetic frustration with increasing 
V gives insight into how a BOW phase might develop, but 
the actual BOW gs for U < U* and V < V C (U) has finite 
P2 that cannot be represented by a HAF with J\, J 2 . 

We consider next the possibility of a physical realiza- 
tion of a BOW phase. Several obstacles were mentioned 
in the Introduction. The material must have quasi-lD 
electronic properties that can be approximated by an ex- 
tended Hubbard model such as Eq. [1] it must be close 
to a CDW instability; and yet it must avoid the Peierls 
instability of half-filled bands. Electrostatic (Madclung) 
energies are inherently 3D. The generalization of cxm in 
Eq. |H to 3D is straightforward and has been applied to 
CT salts [3(| with neutral or ionic gs. The strong 7r- 
acceptor A = TCNQ (tetracyanoquinodimethane) forms 
an extensive series of salts that contain face-to-face stack 
of A~ ion radicals fl3j |. 1:1 salts correspond to a half- 
filled band, and ID Hubbard- type models describe the 
magnetic, optical and electronic properties of TCNQ 
salts. Formally, Eq. [1] refers to electrons in the LUMO 
of TCNQ with i w 0.1 - 0.3 eV. 

We recently pointed out that the A~ A~ stacks in 
Af+TCNQ" salts with M = Na, K, Rb and Cs are 
close to the CDW boundary of AA~ 2 stacks, thereby 
satisfying one condition for a possible BOW phase [3l| . 
Alkali-TCNQ salts are insulators with Coulomb interac- 
tions among localized charges on M + and delocalized 
charges on A~ . Delocalization suggests that the appro- 
priate V m {a) has a > rather than a = for point 
charges, although 3D electrostatic interactions have to be 
taken into account. In any case, we have molecular sites 
in Eq. [1] and the approximation of a rigid lattice with 
uniform spacing must extend to rigid molecules. Alkali- 
TCNQ salts [32| have a phase transition around Td ~ 
300 K to a dimerized stack whose gs is a BOW expected 
in systems with triplet spin excitons [l3| . The A~ stacks 
are regular for T > Td and have inversion symmetry at 
each A~ . A candidate BOW phase then refers to T > Td 
and must have E m ps ksTd in order to justify using the 
gs at K. The molar spin susceptibility, xm(T), of alkali- 
TCNQs is small [13| (w 10" 4 emu) at T d , consistent with 
typical Hubbard- model parameters [3l| . The behavior of 
Xm{T) for T > Td in these salts, however, is not consis- 
tent with a HAF, the early model of choice, or with any 
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Hubbard model with E m = and Xm(0) > 0. Finite E m 
in a stack with equal spacing accounts qualitatively for 
Xm(T) at T > Td- But magnetic data is merely sugges- 
tive of a BOW phase. 

Quasi- ID materials based on stacks of planar molecules 
have strong electron- vibrational coupling to molecular vi- 
brations as well as to lattice modes. The BOW phase is 
an electronic instability that breaks inversion symmetry 
at sites. Lifting inversion symmetry has spectacular vi- 
brational consequences that have motivated Raman and 
infrared studies of dimcrization transitions in many sys- 
tems [33[, including TCNQ salts. The appearance of a 
totally symmetric TCNQ-vibration, polarized along the 
stack, in the ir spectrum at T > Td is evidence for a BOW 
phase. Such a mode was reported [3J| for a CN stretch in 
K-TCNQ and rationalized in terms of dimerization fluc- 
tuations before the suggestion [2j of a BOW phase for 
the EHM. Several coupled modes were reported [35| in 
powder ir spectra at T > Td of several alkali-TCNQs. 
The second form of the rubidium salt, Rb-TCNQ(II), is 
a good candidate for a planned single crystal ir study. It 
has a low (32$ Td « 230 K and a crystal structure (3f| 
with one molecule per unit cell at 300 K. 



In summary, we have characterized the BOW phase of 
half-filled extended Hubbard models with ID potentials 
V m (a) in Eq. [3] The BOW phase (Table 1) for U < U* 
extends from V S (U) where the magnetic gap opens to the 
metallic point V C (U) at the CDW boundary. The kinetic 
energy or gs bond order provides an accurate new esti- 
mate of the critical U* that terminates the BOW phase. 
The SDW/CDW transition for U > U* is a first-order 
transition between two insulating phases. We find that 
the BOW gs has strong charge fluctuations, magnetic 
frustration and a finite two-photon gap. The present 
analysis of competing interactions is limited to electronic 
correlations in a ID model. Both electron-phonon cou- 
pling and electron-molecular- vibration coupling will have 
to be considered to establish that several suggestive ob- 
servations on alkali-TCNQ salts at T > Td can be inter- 
preted as physical realizations of BOW systems. 
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